COMPUTATIONS ON SOME HANKEL MATRICES 



RUIMING ZHANG 

Abstract. In this note, we present the determinant, the inverse and a lower 
bound for the smallest eigenvalue for some Hankel matrices 



1. Introduction 
For each nonnegative integer n, the n-th Hilbert matrix is |4j 

\3 + k + lJ jtk=0 

These matrices are the moment matrices associated the Legendre polynomials. The 
generalized Hilbert matrices, which are also called Hankel matrices, are from the 
generalized moment matrices associated with some more general orthogonal polyno- 
mials. Some interesting questions for Hankel matrices are the determinants, inverses 
and lower bounds for the smallest eigenvalues. In [5] we have developed a general 
method to compute the determinants, inverses and lower bounds for the smallest 
eigenvalues for the generalized moment matrices associated with some orthogonal 
systems (not just limited to orthogonal polynomials). In this note we apply the 
results to some Hankel matrices. The following theorem is adapted from [5] and we 
won't repeat the proof here. 

Theorem 1. Given a probability measure P{dx) on K, for each nonnegative integer 
n, let 

[i n = x n P(dx), 
Jr 

Gn = (/ i i+fc)™ fe=0 

and 

n 

p n {x) = ^a n ^ k x k , n = 0, 1,... 

k=0 

be the orthonormal polynomials, then 

n 

and 

@n = (77,fe)",fe=Q 

2000 Mathematics Subject Classification. Primary 15A09; Secondary 33D45. 
Key words and phrases. Orthogonal Polynomials; Hilbert matrices; Hankel Matrices; Deter- 
minants; Inverse Matrices; Smallest eigenvalue. 



COMPUTATIONS ON SOME HANKEL MATRICES 



2 



with 

n 

£-raax(j,fe) 

Furthermore, if there is a complex number zq with \zo\ — 1 such that for each 
nonnegative integer n, the following sequence 

Q"n,kZ$ > k = 0, 1, . . . 

have the same sign , then the smallest eigenvalue X s of the matrix G n has a lower 
bond 

1 



A, > 



Era=0 

\Pm(z )\ 2 ' 

Remark 2. In the case that all the p m (zo) are real, we could apply the Christoffel- 
Darboux formula to get [TJ [3] 

1n+l,n+l 

A ^ p T" . 

a n ,n {p' n+1 {z )p n (z a ) - p n+1 (z )p' n (z )} 
Recall that the Euler's T(z) is defined as [U El [3] 

T(z) = / x z - 1 e- x dx, > 

Jo 

and it could be analytically extended to a meromorphic function on the complex 
plane. The shifted factorial of z is defined as 

, , T(z + n) 

r(z) 

The hypergeometric function 2F1 is defined as 

v ' n— 

for \z\ < 1. Euler's Beta integral could be evaluated in terms of T(z), 



x 



a-l 



(1 - xf- x dx = , R(a), > 0. 



r(a + /3) 

For any complex number a and < g < 1, we define [TJ [2] 



(a; g)oo = J| a<? m ), (a; <?)m = 



The g-Binomial theorem is 



(az;<l)°c _ y> («;g)fc ,fc | z | < x 



one of its direct consequences is 

00 

0;g)oo = XI 



(-z) fe 



to ^ 
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The confluent g-hypergeometric series i0i 

i / t v {a;g) n q^){-z) n 

i0i (a; o; q, z) = > — r — 

^ (b;q)n(q;q)n 

satisfies the following identity, 

(&/a;g)oo 



0i (a; 6; 3, b/a) 



(M)oo 

Hold |6| < 1 fixed and let a — > oo in the above formula to obtain the Cauchy's 
formula 

oc 

E 



n=0 



{q,b;q) n {b;q) c 
2. Applications 



2.1. Laguerre Polynomials. The Laguerre polynomials {L%(x)}^L are defined 
as [H [2 IE] 



fc=0 

for n > and we assume that 

Lt^x) = 0. 

For any a > — 1, the orthogonal relation for the Laguerre polynomials is 



L«(x)L%(x)^—dx = ^±5 n 
1 (a + 1) ?i! 



for any nonnegative integers to, n where 

1 m = n 



Sn 

Clearly, the n-th moment is 



m n 



J 00 x a+n e~ x dx / 

m ' n = TV i i \ = « + 1 )n, 

r(a + 1) 



and 



e^( x ) = (-i)-J-A— L ^x) 



m 

(aTl)n' 

is the n-th orthonormal polynomial. Apply Theorem [T] with 



_ _/(a+l)„(-n) fe (-iy 

Q>n,k 



n\ {a + l) k k\ 
we get 



det ((a + l)j+fe) < J - fe <„ = J] M a + ' 



fc=0 
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and its inverse matrix is 

,,,,, r i (a + iu-e^j-e), V 

Ua + l^)o<,, fc <„ - ^ £!(a + 1} . (Q + l)kjM ■ 
For a > —1, the smallest eigenvalue of the matrix ((a + l)j + k) Q<j k<n is 



(a + 1)„ 



'0<J,fc<n 

-1 



(n+1)! e^-lJlWf-lj-lW^-l^ft-l)' 
where we applied the formula 



From the the Perron formula [3] 

q sb/2 



for x € C\(0, oo) as n — > oo to get 

- e(-D4 a f(-i) = "°~ le ^ 4 " 1/2) { 



1 + 



as n — > oo, this together with 



T(n + /3) \ \n 



as n — > oo gives 



as ?i — > oo. 



2.2. The Jacobi Polynomials {Pi Q ' W (x)}^ =0 . The Jacobi polynomials (p^'^^a;)} 00 

I J n=0 

are defined as [U El [3] 

™ W " n! 2 X V a+1 '2 

for n > and 

For a, P > — 1 , they satisfy the orthogonal relation 

P^\x)P^\x)w{x)dx = M m „ 

'-1 

for all nonnegative integers n, m where 

w(x) = (l-x^il + x) 13 , 
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and 

2 a+l3+1 T(a + n + 1)T(0 + n + 1) 



K 

Since 



(2n + a + [3 + l)T(a + [3 + n + l)n\ ' 

P ( a ,p) M _ + p ( -n;n + a + (3 + 1 . l + x \ 
n [ ' ~ (-l)"n! 2 1 ^ 13+1 ' 2 j 



we let 



A^\y) = P^\2y-l) 

then, for a, (3 > — 1 we have 



f A^\y)A^\y)w{y)dy = h^S„ 
Jo 



for all nonnegative integers n, m where 

y°(l-^r(a + /3 + 2) 

lyj r(« + i)r(/? + i) ' 

and 

fc( a,/j)_ (" + /3+l)(a + !)„(/? + 1), 



(2n + a + /3 + l)n!(a + /? + 1)„ 
Clearly, the n-th moment is 

(<* + !)„ 



Vn= y n w{y)dy 
Jo 



(a + f3 + 2)„ 
and the n-th orthonormal polynomial is 



(a,/8)/ n = / (2n + a + /? + !)(/? + l) n (a + (3 + l) r 
a " W V (a + /3 + l)(a + l) n n! 

x(-ir 2 F 1 (- n;ri +^+ /3+1 ;y 



Thus, 



(2n + a + (3 + 1)((3 + l)„(a + /3 + 1)„ 



(a + /?+l)(a + l) n n! 
+ a + + l) fc 
X (-l)»(/?+l) fc fcl ' 
and for each nonnegative integer n, we have 



det 



(a + l)j+ fc \ _ -A- (/3 + l) m (a + l) m m! 

> + /? + 2) j+fe J fe <„ ll o (a + (3 + 2) 2m (m + a + (3 + l) m 



COMPUTATIONS ON SOME HANKEL MATRICES 



and 



with 



7j,k 



(a + l) J+k 
(2m + a + (3 + l)(/3 + l) m (a + (3 + 1) T 



E 

m=0 



{a + f3 + l){a + l) m ml 



(—m)j(—m)k(m + a + (3 + l)j(m + a + (3 + l)k 



j\k\{!3 + l)j((3 + l) k 



Then, its smallest eigenvalue of the matrix ^ (a+p+2) k + 

(pi Q,/3) (-3) 



+/3+2) 3 +* J 0<i,fe<n 
2 



has a lower bound 



E 



From the Christoffel-Darboux formula and the formula 



dx 



(n + a + ^+l) (a+ i ;i a + i) 



and 
to get 



£ (^ } (-3)) 



P^(-z) = (-l)^p(^(x) 

+ 1 + n)(a + /3 + 1 +n) 



m=0 



(a + /3 + 1 + 2n)(a + /3 + 2 + 2n)h 
x[(a + /3 + 2 + n )p^+l.«+i) (3)P^' Q 5 (3) 



-(« + /?+! + n)P^(3)pf_r- +1) (3)} 
From the asymptotic formula of Jacobi polynomials to obtain [3] 

Pi-«(3) = {*,(«, ft 3 + 2V2) + O (£) } 



as n — > oo and 



E 

m=0 



(pk^(-3)) ; 
TO) 



0o(a,j9;3 + 2v^)^o(a + l,j9 + l;3 + 2\^) 



r(a+l)r(/3+l)(3 + 2V2) 2 
2r(a + /3 + 2) 



1 + 



as n — > oo. Thus 



( » (P^)(-3)) 2 |" 



>,/3) 



(9 



m=0 "m 



(3 + 2V2) 2 " 



as n — > oo. 
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2.3. The g-Laguerre Polynomials (x; q)}%L with q G (0, 1) and a > — 1. 

The g-Laguerre polynomials {L„ \x; q)}%L are defined as [HE] 

r(a)/ \ _ ^ (q- n ;q) k q^ 1 )q^ k X k 

U [X,Q> (q;q)n ^ (q;q) k (q a+1 ;q) k 

for n > 0, and we assume that 

L^(x;q)=0. 

The moment problem of the g-Laguerre polynomials is indeterminate and one of 
the orthogonality for {L^*\x; q)}^ =0 is 

L^(q k ;q)L^(q k ;q)w ql (q k ;a) = ™ 

, Wi q)n.q 

k— — oo 



for m, n > with 



(9 a+1 ;9)oo(-g;9)oo(-i;g) fc g fe ( Q+1 ) 



(g;<z)oo(-« Q+1 ;'z)oo(-g- Q ;9)oo 

Clearly, the n-th moment is 

oo 
fe=0 

and the orthonormal system is given by 



then, 



u>n,k 



(<z;<?)n (-i)"(<7;9)fe(g Q+1 ;9)fc' 

According to Theorem [IJ the matrix 

(2.i) ((^jffW^ 1 )-^); 

has inverse 



j,k=o 



where 



73,fe fr n (ff;?) m (?;g)i(?;ff)*(ff o+1 ;?)j(9 0,+1 ;?)ik 



m=0 

Its determinant is 



l[(q;q)m(q a+1 ;q)r, 

det ( (a a+1 - a) ^a-( J+ i +i y a(j+k) Y - ^ 

aciy(q ,q) J+k q ) j,k=0 q n(n+l)(4n+6 a +5)/6 

and the smallest eigenvalue has a lower bound is 



-l 



E |4 a) (-i;<z) 



\m=0 
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But 



m[ ,q> V (9; qU i-V m (q;q) k (q a+1 ;q)k ' 



.. { _ ir (<l a+1 ;q)mq m ( m a+1 m+Q+1) 



=(-1)' 



=(-i)' 



f (g a+1 ;g)mg m (g m+a+1 ;g)c 
{q;q) m (g Q+1 ;g)oo 



(g,g a+1 ;g) m ' 

thus the lower bound for the smallest eigenvalue is 

(J? (?,A*) m ) • 

Notice that everything involved here is a rational function of variable q, we may 
replace q by g _1 then apply the relation 

(a;q) n = (-a) n q®(a- 1 ;q- 1 ) n 

and the diagonal matrix 

D n =((-l)^. fe )^ =0 , Dl=I 

to obtain 

n 

det ((g a+1 ;g) i+ fc) <^<„ = a »(«+iwv(«+D(2«+i)/6 n(?;9) fe (g a+1 ;9) fc , 

fc=0 

and 



g ? ' +fc (g a+1 ;g) m (g- m ;g) J -(g- m ;g) fc 
- (<~; g)j (<?; g) fc (g Q+1 ; g),(g Q+1 ; g) fc (g; g) m g( a+1 > 



j,fe=0 

and the smallest eigenvalue has a lower bound 

/ » g 2(?) g m(q + l) \ - 1 

which is bounded again by the absolute constant 



OO • 
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